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A new turbulence model suited for calculating the turbulent Prandtl number as part of 
the solution is presented. The model is based on a set of two equations: one governing the 
variance of the enthalpy and the other governing its dissipation rate. These equations were 
derived from the exact energy equation and thus take into consideration compressibility and 
dissipation terms. The model is used to study two cases involving shock wave/boundary 
layer interaction at Mach 9.22 and Mach 5.0. In general, heat transfer prediction showed 
great improvement over traditional turbulence models where the turbulent Prandtl number 
is assumed constant. It is concluded that using a model that calculates the turbulent Prandtl 
number as part of the solution is the key to bridging the gap between theory and experiment 
for flows dominated by shock wave/boundary layer interactions. 


I. Introduction 

Present simulation of turbulent flows involving shock wave/boundary layer interaction invariably over-estimates 
heat flux by almost a factor of two. 1 One possible reason for such a performance is a result of the fact that the 
turbulence models employed make use of Morkovin’s hypothesis. 2 This hypothesis is valid for non-hypersonic 
Mach numbers and moderate rates of heat transfer. At hypersonic Mach numbers, high rates of heat transfer exist in 
regions where shock wave/boundary layer interactions are important. For such flows, temperature fluctuations, 
which are as important as velocity fluctuations at the higher Mach numbers, play a major role in determining wall 
heat flux and their effects must explicitly be taken into consideration. As a result, one should not expect traditional 
turbulence models to yield accurate results. 

The goal of this investigation is to explore the role of a variable Prandtl number formulation in predicting heat 
flux in flows dominated by strong shock wave/boundary layer interactions. The intended applications involve 
external flows in the absence of combustion such as those encountered in supersonic inlets. This can be achieved by 
adding equations for the temperature (enthalpy) variance and its dissipation rate. Such equations can be derived 
from the exact Navier-Stokes equations. Traditionally, modeled equations (see, for example. Ref. 3,4) are based on 
the low speed energy equation where the pressure gradient term and the term responsible for energy dissipation are 
ignored. It is clear that such assumptions are not valid for hypersonic flows. 

The approach used here is based on the procedure used in deriving the k-'Q model, 5 in which the exact equations 
that governed k, the variance of velocity, and £, the variance of vorticity, were derived and modeled. For the 
variable turbulent Prandtl number, the exact equations that govern the temperature (enthalpy) variance and its 
dissipation rate are derived and modeled term by term. The resulting set of equations are free of damping and wall 
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functions and are coordinate-system independent. Moreover, modeled correlations are tensorially consistent and 
invariant under Galilean transformation. 

Two flat plate experiments are used to determine model constants. The first is the Mach (M) 9.2 experiments of 
Coleman and Stollery 6 , which was conducted in a hypersonic gun tunnel at Imperial College. The second is the 
M=8.3 experiments of Kussoy et al. 7 which were conducted in the Ames 3.5 Foot Hypersonic Wind Tunnel Facility. 
This turned out to be a major undertaking because of the different facilities and instrumentation and because no 
accuracy estimates of heat transfer measurement wave provided. In order to put things in proper perspective, it is 
noted that recent heat transfer measurements 8 on an elliptic cone in the AEDC Tunnel B estimated uncertainties in 
excess of + 10%. The model is validated by recent non-intrusive measurements by Schiilein 9 of flows involving 
shock wave/boundary layer interactions at M = 5. The measurements were carried out at the Ludwig Tube Facility 
at DLR. Oil-film interferometry techniques were used to measure skin friction while an infrared camera system was 
used for heat transfer measurements. 

The recent measurements of Schiilein were a repeat of an earlier experiment 10 which did not include heat transfer 
measurements or skin friction measurements in the separated flow region. Calculations of the earlier experiments 
were carried out by Nance and Hassan 11 using a k-C two-equation and an abbreviated stress model. It was concluded 
in Ref. 10 that there was a need to develop turbulence models capable of predicting the turbulent Prandtl number as 
part of the solution. 

It is shown in this study that a variable Prandtl number formulation results in significant improvement of heat 
transfer predictions in the presence of shock wave/boundary layer interactions. However, the new model has 
insignificant influence on wall pressure and skin friction distributions. 


1. Governing equations 


II. Formulation of the Problem 


The energy equation can be written as 
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and p is the density, h is the enthalpy, p is the pressure, Uj is the velocity and X and (J. are the coefficients of thermal 
conductivity and molecular viscosity. Noting that 


^ = R P ^ + RT^ 
Dt Dt Dt 


( 3 ) 


and using the relation h = C pT and the conservation of mass equation, Eq. (1) can be re-written as 
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p^- = -(y-i)pfc|^- + y + • y =c P ic v 


where C p and C v are the specific heats at constant pressure and constant volume. 
The mean energy equation follows from Eq. (1) as 
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where v is the kinematic viscosity and £ is the enstrophy. 

Equation (4) was the starting point for deriving an equation for the enthalpy variance and its dissipation 
rate. The exact equations are given in the Appendix, while the modeled equations are given as: 
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and where e /, is the rate of dissipation of the enthalpy variance, a is the diffusivity, and 


v, =C U k~ /v^ = C„ Att* 


is the turbulent kinematic viscosity. 
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The modeled equation for the dissipation of enthalpy variance is 
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.. Ch.ii are model constants and are given in Table 1. The turbulent Prandtl number is 

Pr? = v t /a t . 


The choice of a t merits further elaboration. It was indicated in Ref. 4 that experiments in simple shear flows 
showed that the appropriate time scale for temperature fluctuations is proportional to the average of x /, and x k . 
This is the basis for the modeling indicated in Eq. 8. It should be noted that, traditionally 3,4 , the time scale of 
temperature fluctuations is taken as the geometric average of X /, and x k . 


III. Results and Discussion 

The predictions of the model will be compared first with the 15 deg. ramp experiment of Coleman and Strollery. 
In this experiment no flow separation was indicated. The remaining comparisons will be made with Schtilein two- 
dimensional flow measurements for shock angles (3 of 10 and 14 degrees. Flow separation was observed for both of 
these angles. 

The free stream conditions for the experiments of Coleman and Strollery are: M = 9.2, Re = 47 x 10 6 /m, To = 
1070K, Too = 64.5K and T w = 295K.lt has been shown in Ref. 1 1 that use of 241 x 141 Cartesian grid with constant 
spacing in the x(flow) direction and geometric spacing in the y(normal) direction resulted in a grid resolved solution 
and this grid is employee in the present calculations. Figures 1 and 2 compare the pressure distribution for constant 
and variable Prandtl number calculations. As is seen from Fig. 1 the pressure distribution is essentially independent 
of the turbulent Prandtl number. Both constant and variable turbulent Prandtl number calculations overpredict the 
heat flux in the pressure rise region. The variable turbulent Prandtl number calculation is in good agreement with 
experiment in the recovery region while the constant Prandtl number calculations underpredict the heat flux in this 
region. 

The experimental setup of the Schtilein experiment is shown in Fig. 3. A shock generator is mounted on the 
upper wall and the resultant oblique shock wave interacts with the turbulent boundary layer growing on the flat plate 
along the lower wall. The free-stream conditions in the test section were: M = 5, unit Reynolds number = 37 x 
106/m, T 0 = 410K, P 0 = 2.12 MPa and a wall temperature of 300 ± 5K. It was shown in Ref. 11 that a 141 x 141 
Cartesian grid with constant spacing in the x-di recti on and geometric spacings in the v-di recti on resulted in a grid 
resolved solution. A grid of 151 x 141 is used in this investigation. 

The calculations were limited to a region ahead of the point where the reflected shock impinges on the upper 
surface. This assumption makes it possible to use an extrapolation boundary condition at the outflow. Without this 
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assumption, one would be forced to consider the upper wall in the calculations. As is seen below, this assumption 
affected the solution in an adverse manner near the outflow boundary. 

Figures 4-6 compare predictions of surface pressure, wall shear stress and wall heat flux while Figs. 7 and 8 
compare temperature and velocity at x = 576 mm for (3=10 deg. It is seen from Figs. 4 and 5 that the results are 
almost identical for both constant and variable turbulent Prandtl numbers. Both calculations underpredict the 
pressure in the separated region. 

The oil-film interferometry technique cannot be used to determine the extent of the separated region. Instead, 
conventional oil-film visualization was used to deduce the start and end of the separated region. It is seen from Fig. 
5 that the extent of separation is well predicted. However, some discrepancies are noted in predicting the wall shear 
stress in the region. The behavior near the outflow boundary is a result of locating the outflow boundary ahead of 
the point where the reflected shock impinges on the upper boundary. The flow was responding to the compression 
resulting from the reflected shock at the upper wall which is responsible for the increase in wall shear stress. On the 
other hand the calculations were responding to an expanding flow which results in a decrease in shear stress. 

Figure 6 shows that constant Prandtl number calculations overpredict peak heating by a factor of 2. Although 
variable Prandtl number calculations represent a major improvement in estimating the magnitude of the heat flux, 
the choice of the downstream boundary condition appears to have a major influence on heat transfer predictions 
downstream of the separated region. Small unsteadiness was noticed in the calculations but the discrepancy cannot 
be a direct result of that. The fact that heat transfer predictions in the recovery region of the 15 degree ramp were in 
good agreement with experiment suggest that the outflow boundary condition is partially responsible for the 
discrepancy. 

Figures 7 and 8 show that both temperature and velocity profile are in fair to good agreement with experiment. 
The fact that temperature distribution for the variable Prandtl number is in better agreement with experiment is a 
direct result of the fact that the variable Prandtl number formulation results in more realistic heat flux estimates. 

Figures 9-13 compare similar prediction for the (3 = 14 degree case. Similar remarks can be made regarding this 
case. Experimental measurements show more oscillations in the data and this is reflected in the computations. 

IV. Concluding Remarks 

A new approach has been developed for calculating the turbulent Prandtl number as part of the solution. The 
approach is based on a two-equation model for the enthalpy variance and its dissipation rate. All of the correlations 
that appear in the exact equation that govern the enthalpy variance and its dissipation rate are modeled in order to 
ensure the incorporation of relevant physics into the model equations. 

The new formulation is used to study flows characterized by shock wave/boundary layer interactions. In general 
heat flux calculations showed dramatic improvements while surface pressures and wall shear stress were unaffected 
by the variable Prandtl number formulation. 

The discrepancy in the heat flux calculations of the reflected shock experiments of Schtilein when the turbulent 
Prandtl number formulation is employed can be a result of the outflow boundary condition employed, the numerical 
procedure, a modeling issue, or a combination of all of the above. In spite of the discrepancy, the fact remains, 
however, that the key to bridging the gap between theory and experiment in flows dominated by shock 
wave/boundary layer interactions is to use a variable Prandtl number formulation. 
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VII. Appendix 

Exact Equations for Enthalpy Variance and its Dissipation Rate 
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Fig. 1 Computed and measured pressure distribution, 15 deg ramp 
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Fig. 2 Computed and measured heat flux, 15 deg ramp 
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Fig.3 Schematic of shock-wave/boundary layer interaction experiment. 



Fig.4 Computed and measured wall pressure, (3=10 deg 
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Fig. 5 Computed and measured wall shear stress, (3=10 deg 
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Fig. 7 Computed and measured temperature profile at x=376 mm, (3=10 deg 



Fig. 8 Computed and measured velocity profile at x=376 mm, [3=10 deg 
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Fig. 9 Computed and measured wall pressure, (3=14 deg 



Fig. 10 Computed and measured wall shear stress, (1=14 deg 
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Fig.ll Computed and measured wall heat flux, (1=14 deg 



Fig. 12 Computed and measured temperature profile, at x=376 mm, (3=14 deg 
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